arXiv: 1509.04191 vl [astro-ph.SR] 11 Sep 2015 


Circumstellar discs in X/7-ray binaries: first results 
from the Echelle spectrograph 

R. Zamanov^, K. Stoyanov^, J. Marti^ 

^ Institute of Astronomy and National Astronomical Observatory, Bulgarian Academy of 
Sciences, Tsarighradsko Shose 72, BG-1784 Sofia, Bulgaria 
^ Departamento de Fisica (EPSJ), Universidad de Jaen, Campus Las Lagunillas, A3-420, 

23071, Jaen, Spain 

rkz@astro.bas.bg, kstoyanov@astro.bas.bg, jmarti@ujaen.es 
(Submitted on 12.08.2015; Accepted on 08.09.2015) 

Abstract. Here we report our first spectral observations of Be/X-ray and 7 -ray binaries ob¬ 
tained with the new Echelle spectrograph of the National Astronomical Observatory Rozhen. 
For four objects (LSI-|-61°303, 7 Cas, MWC 148, 4U 2206+54), we report the parameters 
and estimate the sizes of their circumstellar discs using different emission lines (Ha, H/3, H 7 , 
Hel and Fell). For MWC 148, we find that the compact object goes deeply through the disc. 
The flank inflections of Ha can be connected with inner ring formed at the periastron pas¬ 
sage or radiation transfer effects. We point out an intriguing similarity between the optical 
emission lines of the 7 -ray binary MWC 148 and the well known Be star 7 Cas. 

Key words: Stars: emission-line. Be - Stars: winds, outflows - X-rays: binaries - Stars: 
individual: LSI+61303, 7 Cas, MWC 148, 4U 2206+54 

1 Introduction 

The Be-X/ 7 -ray binaries (BeX+B) are systems that consist of a compact object 
orbiting an optical companion that is an OBe star. Oe and Be stars are non¬ 
supergiant fast-rotating 0-type and B-type and luminosity class III-V stars 
which, at some point of their lives, have shown spectral lines in emission (Porter 
& Rivinius 2003; Balona 2000; Slettebak 1988). The best studied lines are those 
of hydrogen (Balmer and Paschen series), but the OBe stars may also show 
He and Fe in emission (see Hanuschik 1996, and references therein). They also 
typically display an amount of infrared excess. The origin of the emission lines 
and infrared excesses in BeX+B is attributed to an equatorial disc, fed from 
material expelled from the rapidly rotating OBe star. During periastron, the 
compact object passes close to this disc, and sometimes may even go through 
it causing a major disruption. A large flow of matter is then captured by the 
compact object. In most cases the compact object is a neutron star detected 
as an X-ray pulsar (see Bildsten et al. 1997, and references therein). The con¬ 
version of the kinetic energy of the infalling matter into radiation powers the 
X-rays. Particle acceleration processes associated with the interaction between 
the compact object and the circumstellar disc can also take place, rendering 
the system as a High-Energy (HE) and even Very-High Energy (VHE) source. 

In this paper, we present the current status of a long-term project aimed to 
monitor the variability of emission lines in a selected sample of OBe binaries 
using high resolution spectroscopy. Two of the targets reported here, MWC 
148 and LSI+61°303, belong to the 7 -ray binary family. Only a handful of 
these objects are currently known where the bulk of their luminosity is radiated 
beyond 1 MeV (see e.g. Dubus 2013). The spectral properties of optical emission 
lines, specially in the Ha region, are good tracers of the physical structure 
of the circumstellar discs around the optical companion. It is expected that 
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our observations will enable to establish the expected correlation between the 
spectral line properties and the HE/VHE emission. 

2 Observations and Objects 

The data reported here are from the new EsPeRo spectrograph. EsPeRo is 
hbre-fed Echelle spectrograph attached to the 2.0 m Ritchey-Chretien-Coude 
(RCC) telescope of the National Astronomical Observatory Rozhen (NAO). 
The spectra are reduced in the standard way including bias removal, flat- 
held correction, wavelength calibration and correction for the Earth’s orbital 
motion. Pre-processing of data and parameter measurements are performed 
using various routines provided by the IRAF software package. The wavelength 
calibration is performed using a Thorium-Argon (ThAr) lamp. The light source 
is a 1.5" diameter PHOTRON Hollow Cathode Lamp. To facilitate future 
data reduction, a few Echelle orders with identihed prominent features of the 
ThAr lamp are given in the Appendix. At 6560 A, the spectrograph provides 
a resolution of 0.056 A/pixel, with the Full Width Half Maximum (FWHM) 
of the ThAr lines being 0.22 A. Similarly, at 4800 A, we obtain a spectral 
resolution of 0.041 A/pixel with a 0.16 A FWHM of the calibration lamp lines. 


Table 1. Journal of observations. In the table are given as follows: object 
name, V band magnitude, spectral type, start of the observations, exposure 
time, signal-to-noise ratio at 6600 A, and system orbital period. 


Name 

V 

Spec. 

Date-Obs. 

Exp-time 

S/N 

Porb [fl] 

7 Cas 

2.4 

BO.5 IVpe 

2014-01-14 00:43 

15 min 

90 

203.59 

MWC 148 

9.1 

BO Vpe 

2014-03-13 20:03 

60 min 

65 

321 

LSI+61°303 

10.8 

BO Ve 

2014-02-17 19:23 

60 min 

23 

26.496 

4U 2206+54 

9.9 

09.5 Vep 

2014-01-14 16:12 

60 min 

22 

9.568 


The list of observations is given in Table [H The resulting, normalized spec¬ 
tra are presented in Fig.[T]for different emission lines of interest. The measured 
equivalent widths and the distance between the peaks are summarized in Ta¬ 
ble [21 Here we provide a short description of the four binary stars observed at 
the beginning of 2014. 


7 Cas 

7 Cas is a well known, bright {V ~ 2.3 mag, B0.5IVpe) Be star. It is the 
hrst emission line star discovered by one of the founders of the study of stellar 
spectra, Father Angelo Secchi (Secchi 1866). The orbital period is 203.52 d 
and the orbit seems to be circular (Miroshnichenko, Bjorkman & Krugov 2002, 
Nemravova et al. 2012). Harmanec et al. (2000) reported the mass of the donor 
being ~ 13 Mq and the inclination angle i ~ .^5^. The orbital motion is 
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detectable on spectra and implies that the Be star is the primary component of 
a spectroscopic binary system. The compact object has a mass of about 1 Mq, 
appropriate for a white dwarf or a neutron star, but it could also be a normal 
late-type dwarf (Harmanec et al. 2000; Nemravova et al. 2012, Miroshnichenko 
et al. 2002), or a neutron star in the propeller regime (Shrader et al. 2015), or 
a black hole (low mass, slowly rotating, non-magnetized). Robinson & Smith 
(2000) found that the X-ray flux varied with a period P = 1.12277 d, which 
they interpreted as the rotational period of the mass donor. Here we assume 
Ri = 8.7 ± 0.7 R 0 and v sini= 385 km s“^ (Harmanec 2002) for the radius 
and projected rotational velocity of the donor star, respectively. 

MWC 148 

MWC 148 is a BOVpe star which is the optical counterpart of the TeV (tera- 
electron volt) source HESS J0632+057 (Aharonian et al. 2007). The identifica¬ 
tion of 7 -ray source with the optical Be star is a result of the detection of radio 
(Skilton et al. 2009) and X-ray emission associated with the object (Hinton 
et al. 2009; Falcone et al. 2010) with similar spectral indexes and variability 
as observed in 7 —ray binaries PSR B1259-63, LS 5039 and LSI+61°303 (Kar- 
galtsev et al. 2014). Extended and variable radio emission at AU scales, dis¬ 
playing similar morphology as in these three X-ray binaries, has been reported 
(Moldon, Ribo & Paredes 2011). The absence of significant radial velocity 
shifts in the BO star supports a long period > 100 d (Aragona, McSwain, & De 
Becker 2010). Bongiorno et al. (2011) report strong X-ray flares in Swift/XRT 
data, modulated with 321 ± 5 d, which most probably is the binary period. 
Casares et al. (2012) estimated the radius of the primary Ri = 6 — 1C) R© and 
V sinz= 373 km s“^. 

LSI+61°303 

LSI-|-61^303 consists of a massive BOVe star and a compact object. The orbital 
period is 26.4917 ± 0.0025 d, obtained with Bayesian analysis of radio data 
(Gregory, Peracaula & Taylor 1999; Gregory 2002). The nature of the compact 
object remains a mystery even after 37 years of observations over a wide range 
of wavelengths. Most probably it is a neutron star, but it might be a magnetized 
black hole (Punsly 1999) acting as a precessing microblazar (Massi, Ros & 
Zimmermann 2012). Here we assume a primary radius value of Ri = 6.7 R© 
(Grundstrom et al. 2007) and v sini= 349 ± 6 km s“^ (Zamanov et al. 2013). 

4U 2206+54 

4U 2206+54 is a persistent high-mass X-ray binary star at a distance of 2.6 
kpc (Blay et al. 2006). It was first detected as an X-ray source by the Uhuru 
satellite (Giacconi et al. 1972). The mass donor is classified as an 09.5Vp star 
with a higher than normal helium abundance, underfilling its Roche lobe and 
losing mass via a slow but dense stellar wind, flowing at 350 km s“^ (Ribo et al. 
2006). However, there are some metallic lines typical for a later-type spectrum 
(Negueruela &: Reig 2001). The compact object is a neutron star (Torrejon et 



4 


Zamanov, Stoyanov, Marti 


al. 2004) with spin period Pgpin = 5554 ± 9 s (Finger et al. 2010). The X-ray 
spectrum of 4U 2206+54 is typical for a neutron star accreting material onto its 
magnetic poles. The orbital period of the system is Porb = 9.5591 ± 0.0007 d 
(Ribo et al. 2006) and the orbit is elliptical with e = 0.30 (Stoyanov et al. 
2014). For this object we adopt the values Ri = 7.3 R© (Ribo et al. 2006) and 
V sini= 315 ± 70 km s“^ (Blay et al. 2006). 


Table 2. Measured emission lines parameters. 


object 

line 

EW 

AV 

^disc / -^1 

Rdisc 



[A] 

[A] 


[R 0 I 

7 Gas 

Ha 

25.8 

1.3? 

167? 1400? 


H/3 

3.77 

2.16 

33 

290 


H7 

0.86 

2.60 

18.4 

160 


HeIA6678 

0.23 

6.55 

6.8 

60 


HeIA5876 

0.40 

4.25 

12.5 

110 


FeHA6317 

0.39 

4.47 

13.1 

114 


FeHA5316 

0.28 

3.24 

17.7 

154 


FeHA5169 

0.62 

3.26 

16.5 

144 

MWC 148 

Ha 

34.4 

2.46 

44.0 

352 


H/3 

4.44 

2.38 

25.8 

206 


H7 

1.40 

2.60 

17.2 

137 


HeIA6678 

0.27 

6.24 

7.1 

56 


HeIA5876 

0.39 

4.67 

9.8 

78 


FeHA6317 

0.48 

4.43 

12.6 

101 


FeHA5316 

0.55 

4.25 

9.7 

77 


FeHA5169 

1.10 

3.63 

12.5 

100 

LSI+61°303 Ha 

8.47 

6.56 

5.42 

36.3 


H/3 

0.67 

6.62 

2.92 

19.5 

4U 2206+54 Ha 

0.23 

10.61 

1.69 

12.3 


3 Relation between FWHM, EW and v sin i 

Hanuschik (1989) gives the relation for Pfa emission line: 

log[FWHM(Ha)/(1.23 u sini + 70)] =—0.08 log EW + 0.14, (1) 

where both the FWHM and v sini are in km s“^, EW is in A (see also 
Reid & Parker 2012). Eor LSI+61°303 using Coude spectra with S/N-ratio 
> 60, we measure EWHM(Ha) in the range 12.1 — 13.2 A, and estimate 
V sini= 349 ± 6 km s“^. This is in agreement with the previous measure¬ 
ment V sinz=360 km s“^ by Hutchings & Crampton (1981). Eor MWC 148, 
we measure FWHM(Ha) in the range 10.2 — 10.6 A, and estimate v sini 
= 265 ± 5 km s“^, a value considerably lower than v sini = 373 km s“^ 
(Casares et al. 2012). 
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Fig. 1. Emission line profiles (Ha, H/?, Hel, Fell) of four high-mass X-ray 
binaries (7 Cas, MWC 148, 4U 2206+54, and LSI+61°303). Each spectrum 
is normalized to the local continuum. It is remarkable that both 7 Cas and 
MWC 148 show wine-bottle profile in their respective Ha lines. A clear simi¬ 
larity between the emission lines of these two objects is also visible. The Hell 
4686 A line is not detected in any of the four stars. 
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Fig. 2. The sizes of the disc visible in different emission lines in MWC 148 
seen both face-on (left) and edge-on (right). 



Fig. 3. Comparison of the size of the circumstellar disc and the orbit of the 
compact object in the 7 -ray emitting binary MWC 148. 
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4 Circumstellar disc size 


For rotationally dominated profiles, the peak separation can be regarded as a 
measure of the outer radius (Rdisc) of the Ha emitting disc (Huang, 1972), 


\ ^ /Rdisc 
2v s'mi ) \ Ri ) 


( 2 ) 


Eq. [2] relies on the assumptions that: (1) the Be star is rotating critically, 
(2) the circumstellar disc is Keplerian, and (3) that the line profile shape is 
dominated by kinematics, and radiative transfer does not play a role. Rdisc can 
also be calculated as 


Rdi, 


GMi sin^ i 
(0.5 Z\H)2 ’ 


( 3 ) 


where Mi is the mass of the primary, Ri is its radius, and i is the angle 
of inclination to the line of sight. If the parameters (Mi,i?i,sini) are well 
measured Eqs. [2] and [3] should give similar results. The estimated values of 
disk radius {Rdisc/Ri and Rdisc/R q) are given in Table [2j 

The peak separation in different lines is (Hanuschik, Kozok & Kaiser 1988): 


Z\K(H/?) Ri 1.8Z\K(Ha) (4) 

Z\K(H 7 ) Ri 1.2Z\K(H/?) Ri 2.2Z\K(Ha) (5) 

Z\K(EeH) Ri 2Z\K(Ha) ( 6 ) 

Z\K(EeH) Ri l.lZ\K(H/3) (7) 


where Eq. [7] is derived from Eqs. 0] and [5l 

For MWC 148, we calculate the ratio AV (Yell)/ AV {Ha) Ri 1.96 ± 0.15, 
a value similar to the average value for the Be stars. For MWC 148 and 
LSI+61°303, we also calculate AV(Yifi)/AV(Yia) Ri 1.35, a value considerably 
below the average value for the Be stars. 

For 7 Cas and MWC 148, we calculate Z\K(FeH)/Z\K(H/3) Ri 1.5, a value 
greater than the average value for the Be stars. It seems that the relative size 
of their discs differs from those of the classical Be stars, probably due to tidal 
effects caused by the secondary. 

For the X-ray emission of 7 Cas two mechanisms have been put forward, ac¬ 
cretion of circumstellar disc matter onto an orbiting white dwarf, or magnetic 
field interaction between the star and the circumstellar disc (Smith & Robin¬ 
son 1999). Recently, Motch, Lopes de Oliveira, & Smith (2015) proposed that 
in 7 Cas the magnetic field emerges from equatorially condensed subsurface 
convecting layers the thickness of which steeply increases with rotation rate. 
These authors also propose that 7 Cas and its analogs are the most massive 
and closest to critical rotation Be stars. 


5 "Wine bottle "-type emission line 

In Fig.[Tl it is visible that our observations are with sufficiently high-resolution 
and high S/N ratio to detect "wine bottle"-type flank inflections in Ha emis¬ 
sion of 7 Cas and MWC 148. There are two possible interpretation of the so 
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called "wine bottle" emission-line profiles in Be stars: (i) two component disc 
model (Kogure 1969), where flank inflections are interpreted as due to super¬ 
position of two concentric disc components separated by a gap; (ii) due to line 
radiation transfer effects in an homogeneous, single-component, vertically ex¬ 
tended Keplerian disc as suggested by three-dimensional radiative line-transfer 
calculations (Hummel & Dachs 1992). Symmetric flank inflections of the line 
profiles are then caused by interaction of non-coherent scattering of line pho¬ 
tons in direction perpendicular to the disc plane with rotational broadening. 

At this stage, we consider that both are possible reasons for the "wine 
bottle"-type Ha emission profile of MWC 148. 

Shrader et al. (2015) pointed out a curious resemblance between the mean 
20-60 keV X-ray luminosity of 7 Cas and LSI+61°303. Here we point a clear 
similarity between the optical emission lines of 7 Cas and MWC 148 (see 
Fig.IH). Practically all the emission lines of hydrogen, helium. Fell have nearly 
the same relative intensities and profiles. This indicates that the circumstellar 
disc in these two stars are very similar - close size, density, mass, temperature, 
similar inclination to the line of sight, etc. These coincidences render 7 Cas a 
promising target for future TeV observations, and perhaps even a test bench 
for the future understanding of gamma-ray binaries. 

Let us assume that "wine bottle "-type flank inflections in Ha emission of 
MWC 148 are due to an inner disc/ring. Under this hypothesis, we estimate its 
radius as Rring ~ 68 Rq which is close to the distance between the components 
at periastron passage Vper ~ 86 Rq (Casares et al. 2012). 


6 The 7 -ray binary MWC 148: disk, inclination, resonances 

The measured sizes of the circumstellar disc visible in different emission lines 
in MWC 148 are graphically illustrated in Figs. [2] and |3] from edge-on and 
face-on perspectives, respectively. Our spectrum is obtained at orbital phase 
0.84, which is close to the phase of the periastron (ftp^r ~ 0.97 (Casares et al. 
2012 ). 

A part of the BeXyB have elliptical orbits and there are evidences that 
some systems have a misalignment between the spin axis of the star and the 
spin axis of the binary orbit. The eccentricities in these systems are caused by 
a kick to the neutron star during the supernova event that formed the system. 
Such kicks would also give rise to misalignments (Martin, Tout & Pringel 2009). 
The inclination of the circumstellar disc of MWC 148 to the line of sight is 
probably similar to that of 7 Cas {i ~ 45^, Harmanec 2002). Assuming Mi ~ 
16 Mq, M 2 ~ 4 Mq and aism i ~ 77.6 Rq (see Casares et al. 2012), we 
estimate the inclination of the orbit i ~ 47®. It means that in the case of 
MWC 148, the spin axis of the star approximately coincides with the spin axis 
of the binary orbit. 

For the Be/X-ray binaries, the measurements of disc emission fluxes and 
disk radii show a tendency to have a preferred value (and in the most of 
cases to be around this value, for example see the histogram given in Fig. 6 
by Zamanov et al. 2013), in other words to cluster at specified levels. Haigh, 
Coe & Fabregat (2004) and Coe et al. (2006) discussed that this tendency 
is related to the presence of resonances between the disc gas and neutron 
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star orbital periods. These resonances tend to truncate the disc at specific disc 
radii (Okazaki & Negueruela 2001). The corresponding values can be calculated 
using (e.g. Grundstrom et al. 2007): 



(GMi)V^ 

2 7 m 


Porb-) 


( 8 ) 


where G is the gravitational constant and n is the integer number of disc 
gas rotational periods per one orbital period. For the system MWC 148, we 
calculated the following resonance radii: Rs = 124 Rq, Rj = 135.6 Rq, Rq = 
150.3 R©, R 5 = 169.7 R©, R 4 = 196.9 R©, R 3 = 238.5 R©, R 2 = 312.6 R©, 
and Ri = 496.2 R©. Using e = 0.83, Casares et al. (2012) have calculated 
the distance between the components in periastron and apastron to be Tp^r ~ 
86 R© and lapp ~ 935 R© respectively. Following the disc sizes, given in TableHl 
it seems that the Ha-emitting disc is truncated by 2:1 resonance, H/l-emitting 
disc - by 4:1 and Hy-emitting disc - by 8:1 resonance. 

Recently, Moritani et al. (2015) have found Ha, H/3, and Hy line profiles to 
exhibit remarkable short-term variability after the MWC 148 apastron (phases 
0.6-0.7), whereas they show little variation near the periastron. In contrast to 
the Balmer lines, no profile variability is seen in any Fell emission line. In Fig. 
[3] we also plot the orbit of the compact object in addition to the circumstellar 
disc. Here it is visible how, at periastron passage, the compact object deepens 
inside Ha, H/3, Hy-disc emitting material, but it does not enter the disc domain 
where the Fell emission originates. Consequently, we naturally expect that the 
orbital motion of the compact object will cause profile variations in Ha, H/3, 
and Hy, but almost no variations in Fell lines. This behaviour is partially in 
agreement with the results of Moritani et al. (2015). Future investigations are 
required to answer whether Ha, H/3 and Hy also vary at periastron. 

Conclusions: We report our first observations of Be/X-ray and Be/y-ray 
binaries with the new EsPeRo spectrograph, recently commissioned on the 
2.0m RCC telescope at NAO Rozhen, Bulgaria. Our main results at present 
can be summarized as: 

1. We find a curious similarity between the optical emission lines of the well 
known Be star y Cas and y-ray binary MWC 148, which means that the 
circumstellar discs in these two objects have similar physical parameters. 

2. For four objects (LSI+61°303, y Cas, MWC 148, 4U 2206+54), we measure 
emission line parameters and estimate the size of the circumstellar discs. 

3. For MWC 148, our estimates show that at the periastron passage the com¬ 
pact objects deeply penetrates in the circumstellar disc, that can cause 
changes in the Balmer emission line profiles. 

For the purpose of data reduction, we present the identification of a few orders 
of the comparison spectrum of the Thorium-Argon lamp (see Appendix). 
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Ciencia y Empleo of Junta de Andalucia, and FEDER funds. The manuscript was presented 
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Appendix: Thorium-Argon spectrum for wavelength calibration. 
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Fig. 4. Thorium-Argon Spectrum in the vicinity of Ha line (A6562 A) 
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Fig. 5. Thorium-Argon Spectrum around A6500 A (order ^Ha+1) with some 
prominent features identified. 
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Fig. 6. Thorium-Argon Spectrum around A6250 A, order ^Ha -|- 4. 
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Fig. 8. Thorium-Argon Spectrum in the vicinity of H/3 line (A4861 A) 












